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PROPAGATION  IN  SHIELDS 


A.  Effect  of  Crustal  Structure 

A  basic  understanding  of  the  mechanisms  of  propagation  and  the  genera¬ 
tion  of  regional  phases  is  necessary  for  improving  location,  yield  estima¬ 
tion  and  discrimination  using  regional  phases.  Due  to  the  multi-modal 
structure  of  such  signals,  mode  theory  in  layered  media  needs  to  be  used. 


Ever  since  the  identification  by  Press  and  Ewing  (1952)  of  the  short- 

period  surface  wave  L  as  a  distinct  seismic  phase  (or,  more  properly,  a 

8 

sequence  of  phases),  an  effort  has  been  concentrated  on  understanding  its 
origin  and  the  effects  of  earth  structure  on  its  propagation  and  attenuation. 
In  particular,  it  was  necessary  to  explain  why  was  apparently  unable  to 
travel  through  even  small  segments  of  oceanic  lithosphere  and  to  determine 
whether,  as  it  was  sometimes  conjectured,  was  a  channel  wave  whose  occur¬ 
rence  indicated  the  existence  of  a  high-velocity  lid  trapping  energy  in  a  low- 
velocity  zone  (LVZ) .  A  major  step  towards  understanding  the  nature  of  L 

8 

was  taken  by  Oliver  and  Ewing  (1957),  who  noted  that  the  group  velocities 

of  prominent  maxima  in  the  L  wavetrain  correspond  to  those  of  higher-mode 

8 

Love  waves.  This  identification  of  L  as  a  superposition  of  short-period 

8 

higher-mode  Love  waves  was  reinforced  convincingly  by  Knopoff  et  al  (1973), 
who  demonstrated  by  means  of  synthetic  seismograms  that  realistic  seismic 
source  mechanisms  and  earth  structures  could  produce  higher-mode  Love  waves 
of  sufficient  amplitude  at  the  proper  frequencies  to  account  for  the  obser¬ 
vations  of  L  ,  at  least  to  first  order, 
g 

We  have  undertaken  to  expand  upon  the  modal  superposition  technique  in 
order  to  answer  some  questions  about  L^  propagation  which  ought  to  be  under¬ 
stood  if  L  is  to  be  used  in  a  program  of  seismic  discrimination  at  regional 
8 

distances.  In  particular,  it  is  well  known  that  L  not  only  fails  to  propa- 

8 

gate  through  the  ocean  but  also  attenuates  rapidly  in  certain  continental 

regions.  It  is  also  known  that  the  amplitude  of  L  ,  as  well  as  its  group 

8 

velocity,  is  strongly  Influenced  by  surflcial  geology.  In  order  to  interpret 

observations  of  L  in  terms  of  the  source  mechanism,  as  well  as  to  choose  the 
8 

best  sites  for  locating  future  seismic  stations,  it  is  important  to  have  a 
quantitative  estimate  of  these  path  effects.  We  should  also  like  to  know, 
for  a  given  source-to-receiver  earth  model,  how  large  a  variability  in  the 
characteristics  of  L^  can  be  anticipated  from  earthquakes  of  different  source 
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mechanisms  and  of  different  depths.  It  is  especially  important  to  know  whether 
certain  earthquakes  have  some  L  signal  characteristics  which  can  be  used  as 
diagnostics  in  discrimination  studies.  Our  study  of  synthetic  seismograms 
has  been  aimed  at  investigating  these  questions. 


In  order  to  produce  the  synthetic  seismograms  by  modal  superposition, 
it  is  necessary  first  to  assume  a  flat,  plane-parallel  earth  model  and  to 
calculate  the  dispersion  relation  (i.e.,  phase  velocity  as  a  function  of 
frequency)  for  every  Love-wave  mode  of  this  model.  The  results  of  a  typical 
calculation  are  shown  in  the  accompanying  Figure  1.  The  high-velocity  cut-off 
is  the  shear-wave  velocity  at  the  Moho,  since  higher  velocities  correspond 
not  to  Lg  but  to  the  mantle  channel  wave  (Stephens  and  Isacks,  1977).  In 
the  Interest  of  computational  efficiency,  synthetics  were  sampled  at  five 
points  per  second,  and  thus  no  frequencies  greater  than  the  folding  frequency  of 
2.5  Hz  were  considered;  this  cut-off  results  in  the  inclusion  of  twenty-six 
Love-wave  modes  for  this  particular  earth  model  (Canadian  shield  with  no 
surficial  sediments) .  It  is  important  to  note  the  strong  curvature  in  the 
graphs  of  the  individual  modes  at  certain  frequencies;  for  earth  models  with 
velocity  contrasts  greater  than  those  in  a  shield  model,  a  curvature  can  be 
very  strong  and,  for  models  with  a  crustal  LVZ,  the  modes  osculate.  Because 
the  group  velocity  curves  are  related  by  differentiation  to  the  phase  velocity 
curves,  it  is  necessary  that  the  phase  velocity  curves  be  approximated  accurately 
in  the  regions  of  strong  curvature  or  else  the  synthetics  will  be  marred  by 
spurious  arrival  times,  a  problem  noted  by  Helmberger  and  Harkrider  (1980).  In 
order  to  lessen  the  effects  of  this  problem  we  have  recently  revised  our  program 
to  utilize  a  four-point  Lagrangian  interpolation  for  the  phase  velocity  at  those 
discrete  Fourier  frequencies  used  in  the  synthesis.  Also  shown  are  the 
particle  displacements  as  a  function  of  depth  for  certain  modes  at  a  parti¬ 
cular  frequency  in  Figure  2.  It  is  seen  that  most  of  the  energy  in  the 
fundamental  mode  is  confined  to  the  top  few  kilometers  of  the  crust,  a 
phenomenon  which  will  be  important  in  considering  the  effect  of  a  surficial 
layer  of  low-velocity  sediments.  For  certain  frequencies  and  certain  modes, 
the  displacement  (and  the  tangential  stress)  eigenfunctions  are  not  well  behaved 
in  the  bottom  layer  of  propagation  above  the  attenuating  half-space.  This 
behavior  is  due  to  computational  inaccuracy  in  the  cumulative  multiplication  of 
layer  matrices,  a  phenomenon  well  known  to  all  investigators  in  this  field. 
Finally,  there  is  shown  an  example  of  the  superposition  of  the  twenty-six 
Love-wave  modes  to  produce  L  in  Figure  3. 
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Figure  2.  Particle  displacement  (normalized  to  the  surface 
displacement)  as  a  function  of  depth  for  seven 
selected  Love-wave  modes  at  a  particular 
frequency. 
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Figure  3.  Dispersion  relations  for  26  Love-wave  modes 


r 


Synthetics  have  been  produced  for  a  variety  of  double-couple  source 
mechanisms  located  at  various  depths  and  observed  at  various  distances.  An 
especially  significant  case  which  has  been  treated  is  that  of  a  linear  array 
of  closely  spaced  instruments.  Such  an  array  is  capable  of  measuring  L 

8 

phase  velocity  by  means  of  frequency-wavenumber  processing  (Mrazek  et  al, 

1980) .  Preliminary  results  indicate  that  the  phase  velocity  is  greater  (at 
a  given  frequency)  for  deeper  events,  since  they  preferentially  excite  the 
higher,  faster  modes.  If  a  more  detailed  investigation  which  is  currently 
underway  upholds  this  result,  then  would  in  fact  be  shown  to  be  a  useful 
discriminant  for  earthquakes  occurring  at  depths  of  greater  than  a  few  kilo¬ 
meters.  Also  being  examined  is  the  decay  of  amplitude  with  distance  for  path 
lengths  of  up  to  1800  km.  It  is  important  to  determine  whether  this  decay 
is  sensitive  to  the  source  mechanism,  since  the  surface-wave  excitation  is 
highly  variable,  depending  upon  the  spatial  orientation  of  the  double-couple. 
Not  only  the  amplitude  of  the  L  wavetrain  envelope  but  also  its  temporal 
extent  is  being  compared,  at  certain  distances,  with  observations.  It  appears 
that  the  modal  superposition  predicts  a  wavetrain  of  shorter  duration  than  is 

observed.  The  late-arriving  energy  in  the  observed  L  envelopes  is  probably 

8 

due  to  waves  which  have  been  scattered  by  topographic  and  geologic  inhomo¬ 
geneities  (on  the  scale  of  an  wavelength)  and  which  have  undergone  conver¬ 
sion  to  (slower)  Rayleigh  modes.  Energy  which  is  transmitted  in  the  form  of 
Rayleigh  modes  not  only  lengthens  the  coda  of  L  but  also  gives  rise  to 
particle  displacements  in  the  vertical  and  radial  directions.  This  phenomenon 
is  currently  not  modeled  by  the  synthetics. 

The  Love-wave  dispersion  and  eigenfunction  calculations  have  been  com¬ 
pleted  for  a  Canadian  shield  model  topped  by  a  2-km  thick  sedimentary  layer. 
Synthetics  will  be  run  for  this  second  model  when  analysis  is  completed  of 
those  which  have  been  produced  using  the  hard  rock  model.  Another  model  is 
being  prepared  for  which  the  anelasticity  Q  is  depth-dependent  within  the 
crust.  This  will  result  in  a  surface-wave  attenuation  which  is  both  frequency- 
and  mode-dependent,  and  it  will  be  important  to  determine  what  effect  this  has 
upon  the  amplitude-distance  relation  and  upon  the  F-K  plots  which  would  be 
measured  at  various  distances.  We  note  that  Knopoff  et  al  (1974),  examining 
the  mantle  LVZ  by  means  of  synthesis  at  longer  wave-lengths  than  we  are 
considering,  invokes  a  depth-dependent  Q  as  the  explanation  of  why  L  fails  to 
propagate  in  an  oceanic  earth  model. 
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In  order  to  account  for  the  energy  which  arrives  late  in  the  coda 
in  the  form  of  Rayleigh  modes,  work  is  underway  to  develop  a  Rayleigh-wave 
counterpart  for  the  Love-wave  modal  superposition.  To  this  end,  a  Rayleigh- 
wave  dispersion  program  which  utilizes  the  algorithm  of  Abo-Zena  (1979)  which 
is  more  accurate  for  high  frequencies  and  higher  modes  than  are  the  conventional 
algorithms  has  been  programmed.  The  numerical  difficulties  arising  in  the 
computation  of  the  kinetic  energy  integral  for  the  layered  halfspace  and  in 
the  computation,  via  propagator  matrices,  of  the  vertical  and  radial  dis¬ 
placement  and  stress  eigenfunctions  remain  to  be  addressed. 

D.  W.  Rivers 
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B.  Site  Effects  on  Regional  Phases 

Since  yield  estimation  of  regional  events  Is  based  on  wave  amplitudes, 
and  since  envelope  shapes  may  be  diagnostic  of  source  depth,  the  effect  of 
geological  structures  under  the  observing  sites  on  these  signal  properties 
must  be  evaluated. 


In  our  previous  work  at  NTS  we  have  shown  that  at  sites  underlain  by  low 
velocity  materials,  the  amplitudes  of  the  regional  phases  and  were  ampli¬ 
fied  by  factors  up  to  10,  and  the  envelope  shapes  of  both  phases  also  changed, 
resulting  In  a  prolongation  of  wavetralns  at  those  sites  (Barker  et  al,  1979). 

The  observed  site  effect  for  L  and  P  appear,  therefore,  to  be  several  times 

8  8 

greater  than  that  for  teleselsmlc  P  and  would  seriously  affect  any  yield 
estimates  based  on  measurements  on  these  phases.  The  sites  at  NTS  were 
closely  spaced  such  that  differences  in  paths,  as  manifested  in  multlpathing, 
source  radiation  patterns  and  varying  geophysical  properties  along  the  paths, 
do  not  seriously  affect  the  conclusions  about  site  effects.  To  further  inves¬ 
tigate  site  effects,  we  have  undertaken  the  analysis  of  regional  phase  ampli¬ 
tudes  at  LRSM  stations  located  In  the  north-central  region  of  the  US  and  some 
adjoining  regions  of  Canada.  The  stations  analyzed  thus  far  are  shown  in 
the  enclosed  map  (Figure  4)  showing  the  events  used.  The  event  parameters 
are  listed  in  Table  I.  Amplitudes  of  all  Identifiable  regional  phases  were 
read  at  each  station.  In  addition  for  L^,  the  envelope  shape  was  characterized 
by  noting  the  time  of  the  onset  of  the  phase,  the  time  at  the  maximum  and  finally, 
where  the  envelope  decreased  to  one-half  of  the  maximum.  The  amplitudes  used 
In  the  estimation  of  site  effects  were  at  the  maximum  of  the  magnitude. 


-2 


The  L  amplitudes  were  corrected  to  a  conmion  distance  of  300  km  using  the 


8 


A  “  fall-off  rate  appropriate  to  the  EUS.  The  station  sites  were  subdivided  into 
two  groups:  hardrock  sites  and  sedimentary  sites.  This  subdivision  is  also  indicated 
In  Figure  4.  The  distance  corrected  amplitudes  were  also  corrected  for  event 
magnitude  by  dividing  them  with  the  average  of  hardrock  site  amplitudes. 


Histograms  of  these  reduced  amplitudes  of  for  two  types  of  sites  (Figure  5), 


clearly  show  the  effect  of  site  amplification.  The  average  of  the  sedimentary 
site  amplitudes  is  larger  by  a  factor  of  1.34  (about  .13  amplitude  units).  This 
contrast  Is  considerably  less  than  at  NTS.  This  is  explainable  by  the  fact  that 
the  sedimentary  sites  we  examined  do  not  have  the  thick,  extremely  low-velocity 
sediments  or  volcanlcs  present  at  Pahute  Mesa  and  Yucca  Flats,  and  many  sedi¬ 


mentary  sites  are  actually  on  fairly  well  consolidated  rocks.  The  amplitudes 
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STATIONS 


TABLE  1 


Event  List  For  Regional  Phase  Study 


Date 

Origin 

Time 

Location 

Latitude 

Longitude 

Magnitude 

Depth 

15  Jul  62 

11 

59 

21.9 

Montana 

45.  ON 

110. 2W 

25 

23  Jul  62 

06 

05 

18.4 

New  Madrid 

36.  IN 

89. 8W 

18 

22  Oct  62 

05 

03 

03.9 

Montana 

45. 2N 

111.3W 

33 

04  Nov  62 

06 

18 

31.5 

Yellowstone 

44. 3N 

110. 3W 

04  Dec  62 

17 

49 

59.4 

Colorado 

39. 8N 

104, 7W 

33 

28  Dec  62 

10 

01 

23.6 

Montana 

48. 4N 

113. 9W 

33 

30  Jan  63 

05 

05 

59.5 

Montana 

45.  ON 

110.8W 

33 

30  Jan  63 

23 

05 

09.6 

Colorado 

39. 8N 

104. 6W 

33 

16  Feb  63 

03 

01 

41.0 

Montana 

46.  IN 

lll.OW 

33 

25  Feb  63 

18 

45 

16.5 

Wyoming 

42. 6N 

109 . 2W 

33 

08  Mar  63 

00 

14 

15.6 

Q.  Eliza.  Is. 

76. 7N 

94. 9W 

3.9 

33 

20  Jul  63 

00 

11 

35.0 

N .  Yukon 

65. 2N 

133. 7W 

4.6 

03  Aug  63 

00 

37 

50,3 

S.  Illinois 

37.  ON 

88. 8W 

18 

04  Sep  63 

21 

20 

18.5 

Baffin  Is. 

71. 5N 

72. 8W 

4.1 

33 

04  Sep  63 

21 

41 

00.6 

Baffin  Is. 

71. 6N 

73. 5W 

4.4 

33 

06  Sep  63 

01 

46 

13.0 

Baffin  Is. 

71. 5N 

73.  OW 

4.4 

33 

12  Oct  63 

02 

46 

48.1 

Baffin  Is. 

71. 6N 

73.  OW 

4.1 

33 

15  Oct  63 

12 

28 

58.4 

S.  Quebec 

46 . 6N 

77. 6W 

18 

15  Oct  63 

13 

59 

49.7 

S .  Quebec 

46. 3N 

77. 8W 

3.8 

14 

16  Oct  63 

15 

31 

01.8 

S.  New.  Eng. 

42. 5N 

70. 8W 

20 

08  Jan  64 

10 

04 

31.6 

S .  Quebec 

46.  IN 

77. 7W 

3.8 

33 

18  Feb  64 

09 

31 

10.5 

Alabama 

34. 8N 

85. 5W 

15 

18  May  64 

01 

04 

30,5 

N.  Eliza.  Is. 

74. 3N 

97. 4W 

4.1 

15 

27  Aug  64 

09 

53 

51.1 

N.  Yukon 

65. 3N 

133. 8W 

4.6 

33 

21  Oct  64 

07 

38 

31.0 

Hebgen  Lake 

44. 8N 

111.6W 

5.8 

33 

22  Oct  64 

16 

00 

00.0 

SALMON 

33. 4N 

89. 6W 

4.6 

25  Nov  64 

02 

50 

05.0 

W.  Virginia 

37. 4N 

81. 5W 

4.5 

06  Mar  65 

21 

08 

49.9 

Missouri 

37. 4N 

91.  IW 

5.3 

33 

03  Jun  65 

19 

30 

25.7 

Wyoming 

43. 6N 

106. 5W 

4.7 

33 

15  Jul  65 

14 

16 

07.0 

Off  E.  coast 

37. 3N 

74. 4W 

5.1 

19- 


were  similarly  reduced  by  the  regional  fall-off  rate  of  the  phase  that  can  be 

2  A  1 

written  as  (Alsup,  1970)  A  exp(-irfg-^  •  — )  where  the  values  f  =  3Hz  and 
Q  =  1000  were  used.  The  effect  on  P  is  .26  amplitude  units.  For  P  ,  f  =  2.5 

s  s 

was  used  in  the  same  formula;  the  frequencies  used  were  determined  by  measuring 
the  records.  After  correcting  for  event  magnitudes  we  obtained  histograms  for 
the  two  types  of  sites;  these  are  shown  in  Figures  6  and  7  ,  respectively.  The 
site  effect  due  to  sediments  appears  to  be  quite  large,  .33  magnitude  units. 


In  addition  to  effects  on  amplitudes,  the  envelope  shapes  of  are  also 

changed  by  sediments.  The  observed  effect  is  similar  to  that  observed  at  Yucca 

Flats,  that  is,  the  coda  length  of  L  increases  on  sediments  and  the  maximum 

S 

amplitude  is  also  reached  at  a  later  time.  This  is  illustrated  by  Figure  8. 
This  figure  shows  plots  of  the  arrival  times  of  the  envelope  maximum  and  the 
1/2  max  point,  connected  by  lines,  as  functions  of  epicentral  distance.  The 
two  kinds  of  s)nnbols  denote  the  two  types  of  sites.  The  tendency  is  evident 
for  the  sedimentary  sites  to  have  longer  L  codas  as  well  as  delayed  arrivals 
of  the  maximum  of  the  envelope. 


The  effect  on  envelope  shapes  cannot,  at  present,  be  explained  by 
resonances  in  small  local  structures  as  at  NTS  since  most  of  our  sedimentary 
sites  are  located  in  broad  sedimentary  basins. 


Z.  A.  Der 
A.  O'Donnell 
M.  Marshall 
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Figure  8.  Arrival  times  of  the  L  maxima,  the  1/2  max  point 

(connected  by  bars)  at^hard  rock  (circle)  and  sedimentary 
(+  sign)  sites  plotted  against  distance.  The  L  maxima 
arrive  later  and  the  coda  Is  longer  at  sedimentary  sites. 
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C.  Phase  Velocity  Versus  Depth  Studies 


Apparent  phase  velocities  at  the  maximum  of  appear  to  be  diagnostic 
of  the  source  depth,  and  thus  have  a  discrimination  potential. 

In  an  early  study.  Barley  (1978)  showed  that  apparent  phase  velocities 

measured  across  the  Yellowknife  array  tended  to  increase  with  Increasing 

source  depth.  A  ray  theory  explanation  of  this  phenomenon  is  that  deep 

sources  are  associated  with  rays  of  deeper  penetration  in  the  earth,  with 

correspondingly  higher  apparent  velocities  along  the  surface.  Another 

manifestation  of  this  Idea,  the  change  of  envelope  shapes  with  depth, 

deeper  rays  being  associated  with  faster  group  velocities  was  demonstrated 

by  Noponen  and  Burnetti  (1980).  To  test  this  idea,  we  have  measured  apparent 

phase  velocities  of  L  across  the  Alaskan  array  BFAK.  Although  the  amount  of 

8 

data  processed  is  small,  it  tends  to  confirm  Barley's  findings.  In  Figure  9  we 
show  some  phase  velocity  measurements  versus  source  depth  at  the  Alaskan  array 
BFAK.  The  phase  velocities  were  taken  from  F-K  spectra,  with  fairly  high  F  values 
indicating  good  coherence  across  the  array.  Figure  10  shows  some  phase  velocities 
computed  from  synthetic  seismograms  computed  from  various  source  depths,  also 
showing  the  same  tendency.  We  are  in  the  process  of  collecting  a  much  larger 
data  set  to  investigate  the  stability  of  the  method  observationally  as  well 
as  theoretically.  Clearly,  in  order  to  be  usable,  the  sensitivity  of  the 
method  to  various  crustal  structures  and  source  mechanisms  must  be  known. 

Z.  A.  Der 
P.  Klouda 


D.  Scattering  and  Q  Effects 

Envelope  shapes  of  are  likely  depth  discriminants.  Site  effects  or  the 

time  domain  shapes  of  L  are  caused  by  resonances  in  local,  two  and  three 

8 

dimensional  geological  structures  and  by  scattering  at  near  surface  inhomo¬ 
geneities.  The  relative  importance  of  these  effects  must  be  understood. 

Codas  of  near  earthquakes  have  been  described  in  terms  of  backscattering 
from  random  inhomogeneities,  Aki  (1969)  and  Aki  and  Chouet  (1975).  In  our 
previous  work,  we  described  the  site  amplification  and  coda  lengthening  in 
terms  of  resonance  in  a  basin  filled  with  low-velocity  sediments.  Band-pass 
filtering  of  waves  showed  that  the  coda  length  depends  on  frequency. 

Figure  11  shows  an  example  of  this  at  Yucca  Flats.  Since  all  sites  show  this 
phenomenon,  it  appears  that  there  is  a  contribution  to  the  energy  in  the  L 
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PHASE  VELOCITY  (km/sec) 


DEPTH  (km) 


Figure  10.  Apparent  phase  velocities  (c)  of  synthetic  Lg  (transverse) 
seismograms  versus  source  depth  illustrating  the  principle 
of  dependence  of  c  on  source  depth.  This  figure  is  not 
directly  comparable  to  Figure  9  since  the  component  of 
motion  as  well  as  the  applicable  crustal  structure  is 
different. 


wavetralns  from  scattered  waves  that  should  be  better  understood  If  L  Is 

g 

used  for  yield  estimation.  To  evaluate  this  effect,  we  have  developed 
computer  programs  that: 


Simulate  the  resonances  in  a  2D  basin  structure  with  variable 
velocity  and  Q  structure  using  finite  difference  methods. 

Predict  envelopes  in  scattered  wavetralns  with  various  assumptions 
about  Q  structure,  scattered  and  primary  wave  velocities.  The 
approach  is  similar  to  that  applied  by  Aki  (1969). 

Compute  synthetic  seismograms  using  mode  theory.  The  SH 
(Love)  case  is  currently  operational.  Programs  to  compute 
Rayleigh  mode  synthetics  are  being  developed  currently.  Details 
of  this  development  effort  are  given  in  a  separate  section  of  this 
report. 


Observations  of  L  ,  P  and  P  amplitude  anomalies  and  changes  in  wave 
g  g  ^ 

envelopes  will  be  modeled  at  the  shield  sites  being  studied,  and  compared 
to  synthetic,  finite  difference  and  scattering  models. 


Z.  A.  Der 
A.  O’Donnell 
T.  W.  McElfresh 
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E.  Amplitude-Distance  Relationships  for  Regional  Phases  Based  on  Analysis 
of  Soviet  Data 

Amplitude  distance  relationships  of  the  various  regional  phases  must  be 
known  for  yield  estimation  and  discrimination  using  amplitude  ratios  of  such 
phases . 

Perhaps  the  most  useful  single  source  of  Russian  data  on  regional  phases 
is  the  monograph  by  Antonova  et  al  (1978).  Other  sources  used  in  this  study 
are  Nersesov  and  Rautlan  (1964),  Ruzalkln  et  al  (1977),  and  Shishkevish  (1979). 

Most  of  the  data  In  these  studies  comes  from  selsmographlc  stations  In  and  near  the 
Tien  Shan  region  and  in  the  region  east  of  70*E,  bordering  China  and  Mongolia. 

Most  of  this  area  has  generally  been  recognized  as  one  of  efficient  propagation 
of  Lg  (Gupta  et  al,  1980;  Ruzalkln  et  al,  1977). 

The  amplitude-distance  curves  of  all  regional  phases  (Pn,  Pg,  Sn  and  Lg) 
are  highly  variable  from  one  region  to  another  and  even  from  one  direction 
to  another.  Nersesov  and  Rautlan  used  data  from  a  number  of  stations  In 
the  region  between  the  Pamirs  and  the  Lena  River  and  found  the  spatial 
attenuation  of  all  regional  phases  to  vary  significantly  both  as  a  function 
of  location  and  direction  (see  Figures  11  through  14,  Nersesov  and  Rautlan, 

1964).  Amplitude-distance  relationships  for  different  frequency  bands  at  a 
station  In  southeastern  Kazakh,  central  Asia  again  showed  significant 
variations  with  both  azimuthal  direction  and  average  passband  frequency 
(Antonova  et  al,  1978).  The  amplitude  variations  in  all  these  studies  are 
relatively  the  least  for  the  phases  Lg  and  Pg  although  scatter  by  a  factor 
of  about  2  is  rather  common  even  for  these  two  phases.  The  most  stable  and 
the  largest  amplitude  phase  Is  generally  Lg,  observed  to  distances  of  3000  km 
or  greater.  The  next  most  stable  and  large  amplitude  arrival  is  Pg  recorded 
to  a  distance  of  only  about  1000  km. 

Figure  12  shows  composite  or  average  amplitude  curves  of  Lg,  Pg,  Sn  and 
Pn  from  40  crustal  earthquakes  originating  in  North  Tien  Shan  and  recorded  by 
selsmographlc  stations  equipped  with  SKM-3  instruments  in  the  same  region 
(Shishkevish,  1979).  The  SKM-3  seismograph  is  a  short-period  system  with  flat 
response  in  the  period  range  0.1  sec  to  about  1.2  sec  and  peak  response  at 
about  1.4  sec.  The  peak  amplitudes  have  been  normalized  to  a  standard 
magnitude  and  distance.  The  two  largest-amplltude  phases  are  Lg  and  Pg  and, 
on  the  average,  Lg  Is  about  three  times  as  big  as  Pg. 
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.  Average  amplitude  curves  of  regional  phases  Lg,  Pg, 
and  Pjj  for  the  North  Tien  Shan  region  (after  figure 
Shishkevish,  1979).  The  two  thin  curves  represent 
theoretical  least-squares  fit  to  the  observed  data. 
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Figure  13.  Spectral  amplitude-distance  curves  o 
Northeast  direction  (after  Figure  20 
1978) .  The  thin  lines  represent  thei 
squares  fit  to  the  observed  data. 
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Spectral  amplitude-distance  curves  of  Lg  along  the 
Northeast  direction  (after  Figure  23,  Antonova  et  al 
1978) .  The  thin  lines  represent  theoretical  least- 
squares  fit  to  the  observed  data. 


Lg  may  be  regarded  as  a  higher-mode  surface  wave  traveling  with  minimum 
group  velocity  (Nuttll,  1973)  so  that  the  time-domain  amplitude,  A  varies 
with  distance,  A  as 

-i  -i  -ya 

A  =  KA  3  (sin  A)  2  e  ^ 

(1) 

where  A  is  In  degrees,  Y  Is  the  coefficient  of  anelastlc  attenuation  and  K  Is  a 
constant.  Taking  logarithms  to  base  10,  equation  (1)  gives 


log  A  +  f(A)  =  log  K  -  yA  log  e  (2) 

where  f(A)  Is  a  known  function  of  A.  Equation  (2),  a  linear  equation  suitable 
for  least-squares  curve  fitting,  yielded  a  Y  value  of  about  0.27  per  degree 
for  the  Lg  amplitudes  In  Figure  12.  Recent  data  suggest  that  Pg  may  also  be 
regarded  as  a  dispersed  surface  wavetraln  so  that  equation  (1)  may  be  applicable 
(e.g.  Nuttll,  1980).  Using  data  on  Pg  amplitudes  In  Figure  12, a  mean  value  of 
Y  =  0.30  per  degree  was  obtained.  Therefore  the  attenuation  coefficients  for 
Lg  and  Pg,  based  on  data  In  Figure  1,  are  nearly  equal.  These  values  of  y  for 
Lg  lie  somewhere  between  y  =  0.07/deg  for  1-sec-perlod  Lg  In  the  eastern  United 
States  (EUS)  and  y  -  0.6/deg  for  southern  California  (Nuttll,  1973). 

Data  on  amplitude-distance  curves  of  Lg  and  Pg  for  different  frequency 
bands  (Antonova  et  al,  1978)  provide  valuable  Insight  Into  the  generation, 
propagation  and  spectral  characteristics  of  the  two  Important  regional  phases. 
Recordings  of  earthquakes  by  a  seven-channel  ChISS  station  located  at  Talgar 
(approximately  43®N  and  77°E)  In  Tien  Shan  were  used  In  these  studies.  The 
ChISS  station  Is  a  frequency-selection  seismograph  station  In  which  the  vertical 
component  of  ground  velocity  Is  passed  through  a  system  of  band-pass  filters 
so  that  the  amplitude  of  each  recording  Is  proportional  to  the  vertical  com¬ 
ponent  of  ground  velocity  In  a  particular  frequency  band  (Rautlan  et  al,  1978). 
The  eplcentral  locations  of  earthquakes  were  divided  amongst  four  azimuthal 
directions  relative  to  the  Talgar  station:  Northeast  (Zaysan,  Altay,  Sayan, 
Prlbaykal  and  Transbaykal) ,  East  (Dzungaria,  northwest  China  and  Mongolia), 

South  (Pakistan,  Nepal  and  India)  and  West  (South  Tien  Shan,  Pamirs,  Tadzhik, 
Iran,  and  the  Caucasus). 


A  comparison  of  amplitude-distance  curves  of  L  and  P  ,  (Figures  13  and  16 

8  8 

after  Antonova  et  al,  1978)  shows  several  Interesting  features.  For  propagation 


along  the  same  direction,  the  amplitude-distance  curves  of  Pg  and  Lg  have 
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similar  trends.  Figures  13  and  14  show  the  amplitude-distance  curves,  normalized 
to  magnitude  5  earthquakes,  along  the  Northeast  direction  for  three  frequency 
bands  for  Fg  and  Lg,  respectively.  Amplitudes  of  both  phases  attenuate  slowly 
for  eplcentral  distance,  A  less  than  400  km  and  much  more  rapidly  for  A  greater 
than  500  km.  Both  sets  of  curves  show  a  minimum  at  A  about  800  to  900  km. 
Figure  13  and  16  show  results  for  Pg  and  Lg,  respectively  for  earthquakes  along 
the  West  direction.  An  examination  of  Figures  13  to  16  suggests  that  the 
attenuation  of  both  Pg  and  Lg  varies  greatly  with  direction. 

Assuming  the  validity  of  equation  (1),  the  absorption  coefficient,  y 
and  corresponding  values  of  the  quality  factor,  Q  given  by 

Q  =  —  (3) 

where  u  Is  the  group  velocity,  have  been  obtained  for  the  amplitude-distance 
curves  In  Figures  13  to  16.  The  values  of  u  for  Lg  and  Pg  are  taken  to  be 
3.5  km/sec  and  5.7  km/sec,  respectively  (see  Table  III,  Gupta  et  al,  1980). 

The  results  obtained  by  the  least-squares  linear  regression  of  observed  data 
are  given  In  Table  II,  In  which  the  standard  deviation  error  In  y  Is  also 
Indicated.  In  Figures  13  and  14,  data  for  A  <  500  km  Is  rather  erratic  and 
therefore  not  used  In  the  analysis.  Based  on  the  somewhat  limited  data  In 
Table  II,  the  following  Inferences  may  be  drawn:  (1)  Q  Increases  substantially 
with  frequency  for  both  Pg  and  Lg,  for  A  1000  kn,  (2)  Q(Pg)  Is  less  than 
Q(Lg)  by  a  factor  of  2  to  3  for  a  given  frequency  and  direction  of  propagation, 
(3)Q  (Lg)  Increases  substantially  for  A  >  1000  km,  and  (4)  Lg  travels  more 

efficiently  (higher  Q)  along  the  Northeast  direction  than  along  the  West 
direction  whereas  Pg  Is  less  efficient  for  propagation  along  the  Northeast 
direction  than  along  the  West  direction;  this  effect  Is  more  pronounced  for 
higher  frequencies. 

A  strong  Increase  In  Q  with  frequency  has  previously  been  noted  for  Lg 
by  Molnar  et  al  (1976)  for  the  region  between  Talgar  and  Lake  Baykal 
(Shlshkevlsh,  1979).  Using  data  from  underground  nuclear  explosions  at 
Nevada  Testing  Site,  Press  (1964)  reported  Q  values  for  Lg  greater  than  for 
Pg  by  a  factor  of  2  to  3.  Rautlan  et  al  (1978)  pointed  out  that  Soviet 
scientists  have  often  observed  Q  for  S  waves  to  be  greater  than  Q  for  F 
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Figure  15.  Spectral  amplitude-distance  curves  of  Pg  along  the  west 
direction  (after  Figure  20,  Antonova  et  al,  1978).  The 
thin  lines  represent  theoretical  least-squares  fit  to  the 
observed  data. 


A/T  (micron/sec) 


Figure  16.  Spectral  amplitude-distance  curves  of  Lg  along  the  west 
direction  (after  Figure  23,  Antonova  et  al,  1978).  The 
thin  lines  represent  theoretical  least-squares  fit  to  the 
observed  data. 


TABLE  II 


Direction 

NE 


SPATIAL  ATTENUATION  AND  CORRESPONDING  Q  VALUES 


Range  of  eplcentral 

Mean 

Y 

Phase 

distance  (km) 

frequency  (Hz) 

(per  degree) 

Q 

P 

500  to  1000 

0.35 

0.43  ± 

0.15 

50 

g 

0.70 

0.42  ± 

0.18 

100 

1.4 

0.57  + 

0.20 

150 

L 

500  to  1000 

0.35 

0.29  + 

0.09 

120 

g 

0.70 

0.25  ± 

0.12 

280 

1.4 

0.36  ± 

0.13 

390 

L 

500  to  1000 

0.35 

0.04  ± 

0.02 

925 

g 

0.70 

0.14  ± 

0.02 

495 

1.4 

0.26  ± 

0.02 

540 

P 

300  to  1000 

0.35 

0.31  ± 

0.06 

70 

g 

0.70 

0.35  ± 

0.06 

125 

1.4 

0.27  ± 

0.01 

325 

L 

350  to  1000 

0.35 

0.29  ± 

0.07 

120 

g 

0.70 

0.35  ± 

0.04 

200 

1.4 

0.44  ± 

0.06 

320 

L 

350  to  2500 

0.35 

0.06  ± 

0.02 

580 

g 

0.70 

0.13  ± 

0.02 

550 

1.4 

0.19  ± 

0.02 

725 
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waves;  as  this  Is  opposite  to  observations  of  attenuation  In  the  upper  mantle. 
It  Is  a  puzzling  result.  Q(Lg)  Increasing  with  eplcentral  distance  is,  as 
mentioned  earlier,  probably  due  to  deeper  penetration  of  the  crust.  The 
opposite  directional  effects  suggesting  that  an  Increased  efficiency  of  Lg 
may  be  accompanied  by  a  decreased  efficiency  of  Pg  and  vice  versa  Is  very 
hard  to  explain.  However,  the  presence  of  strong  Lg  and  weak  Pg  In  EUS  and 
weak  Lg  and  strong  Pg  In  WUS  are  perhaps  manifestations  of  the  same  puzzling 
feature  of  Pg  and  Lg  phases.  Most  uf  the  observed  data  came  from  a  single 
station  so  that  the  results  should  be  considered  preliminary. 

I.  N.  Gupta 

J.  A.  Burnettl 
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DETECTION  OF  REGIONAL  PHASES 


The  purpose  of  this  study  was  to  devise  a  regional  event  detector 

making  use  of  the  fundamental-mode  surface-wave  processor  developed  by 

Smart  (1977) .  The  study  has  enabled  us  to  explain  the  past  success  of 

the  processor  as  an  azimuth  estimator,  when  applied  to  L  phases. 

g 


A  test  data  set  of  32  events  was  collected  at  RKON,  well  distributed 
azimuthally  about  that  station  (Table  III),  Figure  17.  Signal  and  noise 
spectra  were  computed  for  all  the  phases.  The  noise  spectra  were  found 
not  to  vary  significantly  over  the  seasons.  Figure  18  shows  the  spectra  of 
four  events,  the  two  highest  and  the  two  lowest  peak  frequencies  in  the 
data  set.  It  also  shows  composite  noise  spectra  by  seasons.  It  can  be 
seen  that  the  noise  spectra  peak  below  0.4  Hz,  whereas  all  signal  spectra 
peak  above  that  frequency.  (The  instrument  response  has  not  been  removed 
for  these  spectra.  Also,  below  0.2  Hz  the  spectra  are  probably  system  noise.) 
In  general,  the  low-frequency  signals  were  those  which  transversed  complex 
moimtainous  regions  on  their  way  to  RKON,  e.g.,  events  from  California, 

Oregon,  Washington,  and  British  Columbia. 


A  modified  version  of  the  Smart  processor  was  used  which  estimates 

particle  motion  ellipticity.  Previously  ellipticity  was  an  input  parameter. 

The  frequency  band  examined  was  0.5  and  3.0  Hz.  The  L  window  for  each  case 

8 

was  hand  picked  by  an  analyst.  Good  back  azimuths  were  obtained  in  general: 
the  rms  error  was  6.7°.  Some  of  the  low-frequency  events  yielded  the  least 
accurate  estimates. 


We  attempted  to  find  a  parameter  for  detection,  i.e.,  one  which  separated 

signals  from  the  noise.  The  parameters  investigated  were:  the  F-statistic, 

the  mean  recurrence  period,  and  the  derivative  -z-rz,  where  E  =  error  and  (t  = 

09 

back  azimuth  (see  Smart,  1977).  The  rms  excursions  for  the  model  in  the 
vertical,  radial,  and  transverse  directions,  denoted  by  Z,  R,  and  T,  were  also 
examined. 


The  derivative  did  not  follow  a  discernible  pattern.  The  F-statistic 
and  mean  recurrence  period  were  consistently  low  for  both  noise  and  signal 
and  did  not  separate  them  (see  Figure  19  for  F-statistic  histogram).  The 
best  discriminator  was  the  quantity  T/R,  but  the  separation  of  the  populations 
was  not  clear  enough  to  provide  a  reliable  detection  parameter  (Figure  20). 

It  was  observed  that  T/R  <  2.0,  when  it  occurred,  accompanied  poor  azimuthal 
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17,  Map  showing  the  location  of  RKON  and  of  the  events 


FREQUENCY  OF  OCCURRENCE 


ERROR  IN  AZIMUTH  OBLATENESS 


P  PHASE 

Figure  21.  Oblateness  values  and  error  from  true  azimuth  as  the 
processor  moves  through  successive  128-point  windows, 
for  event  13  (Montana).  Two  bands  are  displayed  here. 
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estimates,  i.e.,  >  10.0“  error,  and  that  the  smaller  T/R,  the  larger  the  azimuthal 
error. 


These  observations  suggested  that  the  figure  traced  by  the  particle  motion 

in  the  horizontal  plane  might  serve  to  Indicate  signal  azimuth.  An  algorithm 

was  coded  to  find  In  each  time  window  the  orientation  In  the  horizontal  plane 

in  which  the  rms  excursion,  in  the  frequency  band  of  interest,  was  the  least. 

The  azimuth  of  that  orientation  was,  on  the  average,  as  good  an  estimate  of 

surface  wave  origin  as  the  back  azimuth  from  the  original  coherent  processor. 

Moreover,  it  was  more  stable.  In  low  signal-to-noise  records,  where  T/R  fell 

below  2.0,  the  estimated  azimuth  was  more  accurate  than  that  of  the  coherent 

processor.  All  this  appears  to  explain  why  the  Smart  processor  has  been 

successful  In  estimating  L  azimuths,  but  unsuccessful  In  detecting  l  p'nases. 

®  S 

The  azimuth  estimations  have  been  controlled  by  the  dominant  transverse  Love 

component,  but  the  F-statistics  have  been  kept  low  by  the  incoherence  between 

the  vertical  and  radial  components,  incoherence  resulting  from  the  mixing  of 

Rayleigh  modes,  the  fundamental  and  higher  modes.  Thus,  the  Smart  processor 

operating  on  waves  functions  as  an  incoherent  processor. 


The  results  suggested  another  use  of  the  elongation  figure  described  by 
the  particle  motion,  that  is,  to  flag  the  arrival  of  P-waves  and  indicate 
their  azimuths.  Time  windows  of  6.4  seconds — an  order  of  magnitude  smaller  than 
those  for  the  waves — were  used  to  search  in  the  vicinity  of  the  expected 
P  arrivals.  The  axis  of  elongation  was  taken  to  indicate  back  azimuth,  and 
the  amount  of  elongation,  that  is,  the  ratio  of  the  long  axis  to  the  short  axis, 
was  taken  as  the  index  of  signal  presence.  Taking  that  ratio  at  2.4  as  the 
detection  threshold,  in  the  band  0.9  to  6.0  Hz,  80%  of  the  signals  were 
detected  with  a  false  alarm  rate  of  28  per  hour.  The  average  azimuthal  error 
was  6.0  degrees,  and  in  the  case  of  the  low-frequency,  low  signal-to-noise  events 
these  P  azimuthal  estimates  were  markedly  more  accurate  than  those  for  the  L 

g 

waves.  Thus,  in  spite  of  the  high  L  /P  ratios,  the  P  wave  has  proved  more 
useful  for  signal  detection  and  for  azimuth  estimation  than  L^,  within  the 
restraints  of  the  present  investigation. 


The  limitations  of  this  study  have  not  permitted  the  conclusion  of  the 

research.  The  optimum  frequency  band,  the  optimum  signal  flag,  and  the  most 

useful  combination  of  P  and  L  information  have  not  been  determined.  Two 

g 

additional  detection  criteria  have  suggested  themselves,  during  this  work. 
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One  is  the  observed  stability  of  orientation,  over  several  time  windows,  of 
the  elongated  particle-motion  envelope,  both  in  the  P-wave  portion  of  the 


record  and  in  the  L  portion.  The  other  is  the  observed  ninety  degree 
S 

rotation  of  the  envelope  as  the  P-wave  passes  and  the  wavetrain  arrives. 
These  items  demand  further  study. 


Finally,  it  appears  from  the  research  carried  out  thus  far,  that  the 

principal  utility  of  particle-motion  processing  lies  not  in  signal  detection, 

where  the  simple  power  detector  is  so  successful,  but  in  azimuth  and  distance 

determination.  Of  the  missed  signals  in  the  P-wave  detection  trials  discussed 

earlier,  the  P-wave  processor  estimated  the  back  azimuth  of  half  of  them  to 

within  5  degrees  of  their  true  back  azimuth,  even  though  it  did  not  "detect" 

them.  One  of  these  signals  was  not  visible  even  to  the  analyst.  Moreover,  the 

elongation  of  the  figure  traced  by  the  reference  particle  has  not  been  fully 

exploited  in  this  study.  The  orientation  of  the  elongation,  measured  in 

3  dimensions  and  not  just  in  the  horizontal  plane,  will  yield  not  only  azimuth 

but  emergence  angle  as  well  which  indicates  the  source  distance,  and  thus 

fixes  location.  The  time  between  P  and  L  ,  measured  from  the  power  rise  at 

S 

the  signal  onset  down  to  the  point  where  the  particle-motion  orientation 
rotates  90  degrees,  also  yields  a  distance  estimate.  Particle  motion  proces¬ 
sing  can  also  be  exploited  simply  for  picking  phases.  All  of  this  potential  must 
be  addressed  in  future  research. 


E.  Smart 
H.  Sproules 


IMPROVED  LOCATION  TECHNIQUES  WITH  REGIONAL  PHASES 


A.  Results  of  Location  Experiments  Using  P^  Arrivals  Using  HYLO  Methods 

Since  regional  phase  travel  time  curves  depend  on  the  details  of 
crust-upper  mantle  structures  it  is  Important  to  evaluate  the  possibility 
of  improving  locations  by  using  existing  knowledge  of  the  crust  and  upper 
mantle. 

The  evaluation  of  location  accuracies  using  P  and  P  arrivals  with 

n  g 

various  regional  and  local  crustal  models  in  the  Western  United  States  (WUS) 
has  shown  that  (Chang  and  Racine,  1980)  the  Herrin  Model  is  an  overall  good 
average  model  for  the  crustal  structure  in  the  WUS.  In  that  study  varying 
crustal  model  for  each  source  receiver  path  resulted  in  only  small  improve¬ 
ments  in  location  errors,  suggesting  that  the  lateral  heterogeneity  in  WUS  is 
rather  strong. 

This  result  suggested  that  an  alternate  approach,  called  the  HYLO  method, 
reported  by  Herrin  and  Taggart  (1962)  might  be  a  better  way  to  allow  for  the 
lateral  heterogeneity.  In  the  HYLO  method,  the  entire  United  States  is  sub¬ 
divided  by  a  grid  of  one  degree  squares  and  an  estimated  P^  velocity  based  on 
the  crustal  refraction  data  (Figure  22)  is  assigned  to  each  square  in  the 
grid.  The  appropriate  P^  travel  time  is  then  determined  by  ray  tracing.  In 
addition,  in  the  present  study  the  appropriate  source  station  model  from 
Chang  and  Racine  (1980)  is  also  assigned  to  each  source-to-station  path  to 
allow  for  variations  in  crustal  thickness  at  source  or  receiver. 

The  location  accuracy  of  the  modified  HYLO  method  was  tested  with  the 

same  12  nuclear  explosions  used  by  Chang  and  Racine.  Location  errors  were 

compared  to  those  obtained  with  the  standard  Herrin  model,  using  P^  only, 

and  P  and  P  . 
n  g 

Comparison  of  location  errors  shown  in  Table  IV  shows  the  HYLO  method 
is  the  worst  of  three  methods.  There  are  several  possible  explanations  for  this 
poor  result.  First,  the  P^  velocities  assigned  from  interpolating  the 
published  contour  lines  may  not  be  adeqviate.  Second,  the  assigned  crustal 
models  may  not  be  adequate.  A  new  try  should  be  made  with  separate  models 
for  source  and  receivers. 


A.  C.  Chang 
H.  Sproules 
J.  A.  Burnetti 
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22.  Estimated  velocity  in  the  United  States,  based  on  data 
from  deep  seismic  soundings,  underground  nuclear  explo¬ 
sions,  and  earthquakes  (compiled  by  E.  Herrin  and  J.  Taggart) 


B.  Results  of  Successive  Locations  Method 

Conventional  methods  for  calculating  seismic  hypocenter  locations 

do  not  adequately  compensate  for  the  Earth's  heterogeneity.  Traditionally, 

the  effect  of  the  heterogeneous  Earth  in  locating  events  has  been  corrected 

by  master  event  techniques.  Last  year,  Chang  et  al  (1980)  introduced  two 

approaches  to  correct  for  structures  at  the  source  and  station  regions, 

and  also  to  balance  the  effect  of  uneven  distribution  of  seismic  stations. 

These  two  methods  resulted  in  improving  location  errors  using  teleseismic 

P  arrivals.  A  similar  approach  to  location  using  regional  phases  by  Chang 

and  Racine  (1980)  showed  that  location  errors  do  not  improve  unless  many 

local  crustal  models  are  utilized.  The  result  reflects  the  fact  that  the 

Earth's  heterogeneity  is  predominantly  in  the  crust  and  upper  mantle.  A 

subsequent  attempt  to  compute  an  average  P  velocity  for  each  source-to- 

n 

station  path,  called  the  HYLO  method  introduced  by  Herrin  and  Taggart  (1962), 
resulted  in  no  improvement  in  location  errors.  Although  the  complex  nature 
of  the  P^  velocity  variations  in  Nevada-Utah-Arizona  area  is  well  known 
(Figure  22  ) ,  our  experiments  showed  that  the  crust  is  far  more  complex  and 
it  is  not  adequate  to  use  crustal  models  or  use  average  P^  velocities  com¬ 
puted  from  the  published  data. 

A  composite  travel  time  chart  (Figure  23  )  of  7  Nevada  Test  Sites  (NTS) 

explosions  in  the  test  data  base  shows  that  both  P  and  P  are  linear;  there- 

fore,  it  is  possible  to  compute  average  P^  and  P^  velocities  for  the  NTS 

area.  However,  subsequent  individual  travel  time  charts  show  P  velocities 

n 

are  slightly  different  from  event  to  event.  Figures  24  to  26 
show  some  typical  charts  for  NTS  events.  This  indicates  that  an  average 
velocity  may  not  work  as  well  as  an  individual  velocity  estimated  for  the 
individual  event,  particularly  in  the  area  as  complex  as  NTS. 

The  successive  locations  method  uses  the  Herrin  crustal  model  to  compute 
a  preliminary  location.  Then  the  program  computes  P^  and  P^  velocities  using 
the  obtained  preliminary  location.  These  estimated  P^  and  P^  will  be  used  to 
locate  the  event  again.  This  successive  approximation  cycle  can  be 
repeated. 

At  least  three  stations  are  needed  to  determine  the  velocity,  but  the 
reliability  of  the  estimated  velocity  depends  strongly  on  the  number  of  inputs. 
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Figure  24.  Travel  time  curves  of  and  Pg  for  PILEDRIVER  event. 
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Figure  25.  Travel  time  curves  of  Pjj  and  Pg  for  DORMOUSE  prime  event. 
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Figure  26.  Travel  time  curves  of  and  Pg  for  PASSAIC  event. 
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Furthermore,  input  stations  with  travel  time  residuals  greater  th£ui  2.0 

seconds  are  excluded  In  the  velocity  determination,  because  a  station  with 

a  large  error  may  significantly  alter  the  estimated  velocity.  Standard  errors 

are  also  estimated  and  they  replace  the  bullt-ln  standard  error  values 

for  each  phase.  Standard  errors  are  used  to  weight  inversely  the  signal 

phase  In  the  location  process.  It  was  found  that  the  ratio  of  P  and  P 

n  g 

standard  errors  are  near  optimum  at  1.0  to  3.0  as  In  the  bullt-ln 
ratio. 

Table  V  shows  the  results  In  terms  of  location  error  vectors  and  origin  time 
errors.  Substantial  Improvements  In  location  errors  can  be  seen.  Origin  time 
errors  and  depth  estimate  errors  are  only  as  good  as  those  using  the  Herrin 
travel  time  tables. 

Romney  et  al  (1962)  demonstrated  that  there  are  clearly  two  distinct  P^ 
velocities  toward  the  East  and  toward  the  Northwest  of  the  GNOME  location. 

The  location  error  of  GNOME  is  7.3  km  with  Herrin  and  7.6  km  with  the  suc¬ 
cessive  location  method.  To  account  for  two  clearly  different  crusts,  an 

option  to  estimate  two  sets  of  P  and  P  velocities  was  created  and  tested. 

g  n 

Using  this  dual  successive  determinations  method  the  GNOME  location  error  was 
reduced  from  7.0  km  to  5.4  km. 


A.  C.  Chang 
J.  A.  Burnettl 
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TABLE  V 

COMPARISON  OF  LOCATION  ERRORS 

SUCCESSIVE  LOCATIONS  METHODS 
DEPTH  FREE 


HERRIN 

SUCCESSIVE  LOCATIONS 

LOC  E^. 

OT.  E^. 

depth!^ 

LOC.  E. 

OT.  E. 

DEPTH 

Faultless 

15.04 

1.4 

0.0 

3.96 

1.0 

3.6 

Rullson 

L.58 

0.0 

0.0 

3.39 

2.6 

0.0 

Passaic 

4.08 

-0.2 

0.0 

2.14 

-0.6 

0.0 

Rockville  Dam 

5.18 

1.5 

0.0 

2.07 

3.9 

13.0 

Dormouse  Prime 

9.09 

0.2 

0.0 

9.87 

0.7 

4.3 

Klickitat 

6.53 

0.1 

0.0 

2.78 

0.7 

8.1 

Bandicoot 

1.24 

0.6 

7.8 

4.58 

0.0 

0.0 

Shoal* 

5.18 

0.0 

0.0 

2.34 

-0.9 

0.0 

Merrimac 

9.58 

0.3 

0.0 

5.16 

0.2 

15.9 

Gasbuggy 

11.43 

0.6 

0.0 

10.04 

2.0 

0.0 

Piledriver 

12.98 

0.0 

0.0 

1.66 

-0.1 

8.9 

Roanoke 

8.24 

1.3 

7.9 

4.20 

1.1 

28.1 

Average 

7.52 

4.35 

*  Depth  restricted 

1.  LOC.  E.:  Error  of  location  vectors  in  km 

2.  OT.  E.;  Error  of  origin  time  estimates,  T  .  -  T  ,  in  sec 

est  true 

3.  when  depth  is  zero,  the  depth  is  restricted  to  0  km 
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C.  Improved  Location  With  Regional  Phases 
Method  of  Simultaneous  Inversions 


The  standard  technique  of  seismic  event  location  Involves  assuming  a 

trial  value  of  the  hypocenter  and  origin  time,  calculating  predicted  arrival 

times  at  a  set  of  stations  for  signals  originating  at  the  trial  hypocenter, 

and  expanding  the  arrival  time  residual  6t*t.  .-t,,^.  according 

observed  calculated 

to 

,  at  .  .  at  .  .  at  .. 

=  ax  a^  ^  ah 


and  then  adding  the  corrections  dT,  dx,  dy,  and  dh  to  the  assumed  values  of 

the  origin  time,  east,  north,  and  depth  coordinates.  In  order  to  evaluate 

this  equation  for  each  residual  6t,  at/ai  Is  set  equal  to  unity  at  all  stations 

and  the  other  three  partial  derivatives  are  calculated  by  differentiating  the 

travel-time  curve  for  the  signal  In  question  (P,  P^,  or  P^)  at  the  proper 

epicenter  to  station  distance  A.  In  order  to  perform  this  differentiation. 

It  Is  necessary  to  know  accurately  the  values  of  the  travel- time.  In  particular, 

for  P  and  P  the  travel-time  curves  are  known  to  be  of  the  form 
n  g 


t  =  T  +  Sp  (h)  +  bp  A 
*^n  n  n 

where  the  coefficients  a  and  b  may  be  calculated  from  the  assumed  values  of 
the  thickness  and  the  P-wave  velocity  for  each  layer  of  the  crust  In  a 
laterally  homogeneous  earth  model.  The  HYLO  method  consists  of  using 
locally  measured,  accurate  values  of  the  thicknesses  and  velocities  for  two- 
layer  models  In  various  regions  of  the  earth  In  order  to  more  accurately 
estimate  the  coefficients  a  and  b  and  hence  the  partial  derivatives  dt/3A 
and  dt/3h.  The  method  of  successive  determinations  assumes  the  thicknesses 
of  the  crustal  layers  to  be  known,  but  It  determines  the  velocities  (and 
hence  the  parameters  b  and  b„  )  by  Iteration.  We  describe  herein  a  third 
method  for  differentiating  the  travel-time  curves,  one  which  assumes  the 
coefficients  a  and  b  to  be  unknown  parameters  which  are  to  be  determined  along 
with  the  hypocentral  coordinates,  T,  x,  y,  and  h.  For  n^  observations  of 
teleselsmlc  P,  n,  observations  of  P  ,  and  n.  observations  of  P  ,  this  method 
results  In  n^^  equations  for  arrival  time  residuals  having  the  form  described 
previously,  n^  equations  of  the  form 
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dT 


j.  . 


j.  51,  . 


It 

3  a, 


‘^®Pr, 


3_t 
3  b, 


and  n.  equations  of  a  similar  form  Involving  a-  and  b  .  The  partial  derl- 

vatlves  In  these  equations  may  be  evaluated  In  a  straightforward  manner  from 

the  linear  equations  for  the  P  and  P  travel  times.  The  travel  times  for 

“8 

the  regional  phases  are  no  longer  explicitly  depth-dependent,  since  all  depth 
dependence  Is  absorbed  within  the  Independent  parameters  apjj  and  ap  .  Without 

O 

teleseismic  data,  it  is  Impossible  to  solve  for  the  depth.  Using  the  determined 

values  of  bpg  and  bp^  however,  one  may  find  the  P-wave  crustal  velocities  and, 

if  one  assumes  the  thickness  to  be  known,  one  may  then  solve  for  the  depth  using 

the  values  found  for  ap  and  ap  .  As  can  be  seen  from  the  linear  equation  for 

n  8 

the  P^  and  P^  travel  times,  it  is  impossible  to  separate  the  origin  time  from  ap^ 
or  ap  ,  so  the  origin  time  can  be  determined  only  if  there  are  teleseismic 
measurements  or  if  there  are  both  P  and  P  measurements. 


n  g 

In  order  to  Implement  the  method  of  simultaneous  inversions,  the  standard 
location  program  LOCATION  has  been  rewritten  to  solve  for  eight,  rather  than  four, 
unknowns.  Preliminary  testing  of  the  revised  program  shows  it  to  be  working 
satisfactorily.  A  detailed  evaluation  of  the  program  is  currently  underway, 
using  the  same  data  base  which  was  used  for  the  method  of  successive  determinations. 


The  coefficients  ap  ,  bp  ,  ap  and  bp  are  the  same  for  all  events  with 
epicenters  close  to  each  other.  They  may  thus  be  determined  more  accurately 
if  data  from  several  nearby  events  are  inverted  simultaneously,  a  process 
known  as  joint  epicenter  determination  (JED) .  The  use  of  n  epicenters  results 
in  a  system  of  equations  in  4n+4  unknowns.  After  the  evaluation  of  the  current 
version  of  the  method  of  simultaneous  inversions  is  completed,  the  method  will 
be  expanded  to  encompass  the  use  of  JED.  The  rewriting  of  program  LOC  was 
carried  out  in  such  a  manner  as  to  facilitate  this  extension,  and  the  program¬ 
ming,  testing,  and  evaluation  of  the  method  is  expected  to  be  completed  by 
mid-June . 


D.  W.  Rivers 
A.  C.  Chang 
J.  A.  Burnettl 
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D.  Location  with  Regional  Data 

The  traditional  determination  of  hypocenters  and  origin  times  of  seismic 
events  Involves  finding  the  best  values  of  the  coordinates  (latitude,  longi¬ 
tude,  depth,  origin  time)  such  that  at  a  network  of  detecting  stations  the 

signal  arrival  time  differences  6t»t.  .-t  ,  ^  ,  are  minimized 

®  observed  calculated 

In  a  least-squares  sense.  In  order  to  carry  out  the  location,  it  is  neces¬ 
sary  to  have  measured  arrival  times  at  four  or  more  stations  In  order  that 
each  of  the  four  coordinates  of  the  hypocenter  can  be  determined  uniquely. 

An  accurate  location,  l.e.,  one  for  which  the  "error  ellipsoid"  formed  by 
the  confidence  limits  surrounding  the  calculated  hypocenter  is  small,  would 
obviously  require  several  more  than  four  measurements  of  arrival  time.  In 
practice,  this  requirement  may  be  difficult  to  meet.  One  case  In  which  It 
may  not  be  met  Is  that  of  a  weak  event  for  which  the  P-wave  signal  cannot  be 
observed  at  teleselsmlc  distances.  Unless  the  area  around  the  epicenter  Is 
densely  populated  by  selsmographlc  stations.  It  Is  unlikely  that  enough 
detections  will  be  reported  by  stations  at  only  regional  distances  to  permit 
an  accurate  determination  to  be  carried  out.  Another  case  In  which  the 
location  may  suffer  from  an  Insufficient  number  of  arrival  time  measurements 
Is  that  rather  common  Instance  In  which,  within  some  given  time  Interval, 
several  stations  In  a  network  have  measured  arrival  times  for  P-waves,  but 
the  signals  are  from  two  or  more  different  events  (Including  local  ones), 
and  It  Is  Impossible  to  associate  unambiguously  each  signal  with  the  appro¬ 
priate  event. 


These  shortcomings  of  conventional  location  may  be  circumvented,  or  at 
least  alleviated,  by  the  use  of  more  Information  than  just  the  arrival  time 
measurements.  Smart  (1977)  has  developed  a  surface-wave  processor  which 
enables  three-component  data  to  be  Inverted  to  yield  the  back  azimuth  of  the 
epicenter  at  the  detecting  station.  The  back-azimuth  measurements  made  by 
this  processor  are  an  excellent  supplement  to  the  arrival-time  data.  Even 
If  no  teleselsmlc  measurements  are  available  for  an  event,  the  surface-wave 
processor  may  still  be  used  since,  for  regional  events,  the  amplitude  of 
Is  comparable  to  that  of  P.  In  the  case  in  which  the  P-wave  detections 
cannot  be  associated  unambiguously  with  a  particular  event,  measurement  of 
the  back  azimuth  will  In  many  cases  resolve  the  ambiguity.  In  particular, 
arrival  time  and  back  azimuth  measurements  at  only  two  stations  would  be 
sufficient  to  determine  the  hypocenter.  Although  this  determination  would 
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be  crude,  it  should  serve  as  an  adequate  trial  value  to  permit  the  proper 
association  of  other  signals  with  it.  Back  azimuth  measurements  are  thus 
useful  for  event  location,  and  we  therefore  have  revised  our  standard 
location  program  LOCATION  to  use  azimuth  as  well  as  arrival  time  data. 

In  the  revised  version  of  LOCATION,  an  initial  trial  value  of  the  hypo- 
center  and  origin  time  is  used  to  calculate  the  predicted  signal  arrival 
time  and  back  azimuth  at  each  station  in  a  given  network.  Residuals  are 
calculated  by  subtracting  the  predicted  values  from  the  observed  values, 
and  these  residuals  are  then  expressed  in  the  form  of  a  Taylor  series  about 
the  trial  coordinates.  For  a  data  base  consisting  of  m  arrival  time  measure¬ 
ments  and  n  back-azimuth  measurements,  this  process  results  in  m  equations 
of  the  form 

fit  =  dT  +  |4  dx  +  l^y  +  l^h 
9x  9y  9h 

and  n  equations  of  the  form 

where  6t  and  6^  are  the  arrival  time  and  back  azimuth  residuals  and  dT,  dx, 
dy,  and  dh  are  the  corrections  which  are  to  be  added  to  the  trial  value  of 
the  origin  time,  east,  north,  and  depth  coordinates.  The  partial  derivatives 
in  the  arrival  time  equation  may  be  calculated  in  the  usual  manner  from  travel- 
time  curves,  and  those  in  the  back  azimuth  equation  may  be  determined  by  a 
straightforward  solution  for  the  appropriate  spherical  triangles.  The 
resulting  system  of  nr+n  equations  may  then  be  solved  for  the  four  unknowns, 
dT,  dx,  dy,  and  dh;  these  values  are  then  added  to  the  assumed  hypocentral 
coordinates,  yielding  an  improved  location.  This  location  is  then  used  as 
a  new  trial  value,  and  the  process  may  be  repeated  through  several  iterations 
until  the  location  converges.  In  practice,  every  residual  6t  ought  to  be 
weighted  by  the  standard  deviation  of  the  appropriate  arrival  time  measure¬ 
ment  (e.g.,  0.2  sec  for  P,  1.0  sec  for  P^,  3.0  sec  for  P^,  etc.)  and  every 
residual  6^  ought  to  be  weighted  by  the  standard  deviation  of  a  back  azimuth 
measurement,  which  Smart  (1977)  finds  to  be  7  degrees.  Since  the  intersection 
of  circles  of  radius  t±1.0  sec  around  each  station  comprises  a  much  smaller 
area  than  does  the  Intersection  of  sectors  of  angle  ^+7  degrees  from  each 
station,  it  is  seen  that  the  use  of  back  azimuth  data  will  offer  very  little 
Improvement  in  the  location  beyond  that  which  would  be  obtained  by  the  use 
of  several  measurements  of  only  the  arrival  time.  Nevertheless,  the  back 
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•slmuth  data  will  be  helpful  when  few  arrival  times  were  measured;  in  fact, 
the  epicenter,  but  not  the  depth  or  origin  time,  may  be  calculated  (crudely) 
from  two  back  azimuth  measurements  alone.  It  may  thus  be  possible  to  use 
this  technique  for  locating  "lonesome  L^"  events. 

The  writing  and  debugging  of  the  revised  LOCATION  program  have  been 
completed.  The  program  was  tested  with  azimuth  inputs  shown  in  Smart  (1977) 
for  SALHON  and  GNOME.  Location  errors  are  417.4  km  for  SALMON  and  202.3  km 
for  GNOME.  This  result  is  approximately  the  same  as  the  location  errors  found 
by  Smart. 

D.  W.  Rivers 
A.  C.  Chang 
J.  A.  Bumettl 
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EVASION 


By  detailed  studies  of  explosions  detonated  In  salt  In  the  USSR 
we  hope  to  better  understand  their  source  mechanism.  This  knowledge 
will  make  It  possible  to  evaluate  possibilities  for  evasion  by  decoupling 
In  the  USSR.  For  example,  If  the  SALMON  RDP  Is  a  satisfactory  fit  to  the 
data  then  one  might  more  reasonably  use  decoupling  values  and  RDF's  from 
STERLING  In  analyses  of  Soviet  evasion. 

One  of  the  principal  tools  which  we  plan  to  use  for  Investigation  of 
the  Soviet  shots  Is  the  maximum  likelihood  estimator  developed  by  Shumway 
and  Blandford  (1977).  To  check  out  Its  application  to  salt  shots  we  have 
applied  It  to  the  Initial  P  data  recorded  from  the  event  SALMON.  Figure  27 
shows  the  likelihood  contoured  as  a  function  of  delay  and  amplitude  of 
the  secondary  arrival.  The  source  spectra  used  was  that  for  SALMON,  see 
e.g.  Blandford  and  Woolson  (1979),  and  the  t*  values  were  computed  from  the 
study  of  Der  and  McElfresh  (1976). 

The  station  plots  are  arranged  In  three  profiles  to  the  Northeast,  North, 
and  Northwest.  We  see  that  at  all  azimuths  a  pP  echo  with  a  delay  of  0.55 
±  0.02  seconds  and  a  reflection  coefficient  of  approximately  -0.4.  The 
confidence  Intervals  on  arrival  time  are  the  population  standard  deviation  of 
times  determined  at  the  16  stations.  Only  at  station  WFMN  does  the  reflection 
appear  to  be  absent. 

We  see  that  there  is  also  a  multipath  arrival  with  amplitude  approxi¬ 
mately  +0.4  and  arrival  times  ranging  from  0.1  to  0.45  seconds  with  a  mean 
and  population  standard  deviation  of  0.38  ±  .08. 

In  addition,  we  occasionally  see  a  negative  amplitude  arrival  at  small 
delays.  This  is  especially  prominent  at  the  four  closest  stations  on  the  North 
profile.  We  interpret  this  as  a  reflection  from  the  interface  between 
the  top  of  the  salt  dome  and  the  overlying  sediments.  This  could  account 
for  Its  erratic  appearance  since  the  salt-sediment  interface  is  doubtless 
more  variable  as  a  function  of  azimuth  and  angle  to  the  vertical  than  is 
the  free  surface.  One  can  be  fairly  sure  that  this  is  not  just  an  artifact 
of  the  data  processing  because  of  the  fact  that  Shumway  and  Blandford  (1977) 
found  no  such  reflections  in  the  analysis  of  nine  other  events.  Such  a 
reflection  would  not  have  been  expectr.d  from  those  events  since  they  have 
no  internal  Interfaces  comparable  in  Impedance  to  the  salt-sediments  inter¬ 
face. 
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Ac  the  top  of  Figure  27  is  a  calculation  performed  on  the  average 
of  all  16  spectra,  averaged  after  the  effects  of  t*  had  been  taken  out 
by  divillng  by  exp(-wt*/2).  We  see  that  Che  main  pP  and  the  multipath 
arrival  can  still  be  seen. 

The  success  of  this  approach  as  applied  to  SALMON  suggests  that  the 
SALMON  RDP  is  appropriate  to  SALMON  as  observed  Celeselsmlcally.  When 
moderate  variations  of  the  SALMON  ROP  are  tried,  pP  cannot  be  detected 
trlth  this  approach.  Thus  Che  success  of  this  method  shows  that  the  SALMON 
ROP  is  appropriate  for  the  event. 

In  the  lower  left-hand  comer  of  Figure  27  we  see  the  same  procedure, 
using  the  SALMON  RDP  and  an  independently  determined  t*,  to  analyze  an 
April  25,  1975  shot  from  Azgir  as  recorded  at  NORSAR. 

R.  R.  Blandford 


SUMMARY 


The  main  findings  of  the  research  reported  can  be  summarized  as 
follows: 


•  The  short-period  surface  wave  has  been  modeled  by  means  of 
synthetic  seismograms  in  order  to  determine  the  effects  on  the  observed 
waveforms  of  propagation-path  structure  and  to  evaluate  the  feasibility 
of  using  as  a  depth  discriminant  for  shallow  events. 

•  Amplitudes  of  regional  phases  P^,  and  are  significantly 
amplified  by  sediments.  This  affects  yield  estimation  from  regional  phases, 
as  well  as  discrimination. 

•  Envelope  shapes  of  L  wavetrains  are  also  changed  by  sediments 

S 

and  this  may  Influence  any  attempt  to  determine  source  depth  from  L  . 

g 

•  Phase  velocities  of  across  small  arrays  are  potentially  useful 
as  depth  discriminants. 

•  Besides  near  surface  sediments,  scattering  and  local  near-surface 
attenuation  are  also  critical  factors  affecting  regional  phase  characteristics. 

•  Analysis  of  Soviet  data  on  amplitude-distance  relationships  of 

regional  phases  has  been  carried  out  to  understand  the  effects  of  frequency 

and  propagation  path.  The  two  larges t-amplitude  phases,  L  and  P  ,  show 

g  g 

strong  directional  effects  and  Q  increasing  with  frequency. 

•  An  attempt  was  made  to  devise  a  regional  event  detector  using  the 

Smart  processor.  Signal  and  noise  spectra  were  analyzed  to  determine 

differences  in  spectral  shape.  Good  back-azimuths  were  obtained  using  the 

coherent  fundamental-mode  processor  on  the  L  phases,  however,  no  reliable 

g 

noise/signal  discrimination  statistic  could  be  found.  A  new  processor  was 

developed,  which  used  only  the  incoherent  Love  component  of  the  signal.  This 

processor  gave  excellent  back-azimuths  for  both  the  L  and  P  waves  of  the  event. 

g 

The  oblateness  of  the  motion,  i.e.,  the  ratio  of  the  major  axis  to  the  minor  axis 
of  the  envelope  of  particle  motion,  is  a  reliable  detector  but  not  necessarily 
better  than  a  power  detector. 

•  The  test  results  of  Herrin  and  Taggart's  HYLO  method  were  about 
9.5  km  In  average  errors.  The  failure  to  Improve  location  accuracies  may 
have  been  caused  by  the  complexity  of  the  WUS  area.  In  the  second  subtask, 
the  successive  location  methods  showed  an  average  location  error  of  about 
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4.3  km.  This  result  was  very  encouraging. 

•  A  technique  has  been  developed  which  performs  event  location 
while  treating  the  slopes  and  Intercepts  of  the  (linear)  P  and  P 
travel  time  curves  as  unknowns,  thereby  enabling  an  epicenter  to  be 
determined  without  reference  to  any  a  priori  local  earth  model. 

•  The  program  modification  to  accept  azimuth  inputs  from 
Smart's  detector  has  been  completed  and  tested.  A  preliminary  test 
shows  the  location  with  azimuths  only  is  capable  of  giving  epicenter 
estimates  which  are  comparable  with  locations  obtained  by  Smart  using  another 
network  method. 

•  A  method  has  been  developed  for  increasing  the  amount  of  data 
available  for  use  in  event  location  by  incorporating  into  the  data  base 
measurements  of  back-azimuth  made  by  Smart’s  (1977)  detector. 

•  The  success  of  the  application  of  the  maximum-likelihood 
determination  techniques  of  Shumway  and  Blandford  (1977)  to  SALMON  and 
to  an  Azgir  event  of  April  25,  1975,  suggests  that  the  latter  events  RDP 
may  be  adequately  modeled  by  the  SALMON  RDP. 
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